Numerical analysis of MHD axisymmetric rotating Bodewadt rheology under viscous dissipation and ohmic heating effects

In present research manuscript, analysis is presented for the influences of heat transition in a bodewadt flow over a penetrable disk numerically. Estimation parameters in current mathematical flow model include magnetic field parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(0.1 \le M \le 1.2),$$\end{document}(0.1≤M≤1.2), wall suction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(1.7 \le A \le 6.7),$$\end{document}(1.7≤A≤6.7), prandtl number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(0.2 \le \Pr \le 5.0),$$\end{document}(0.2≤Pr≤5.0), heat generation/absorption \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${( - }0.9 \le Q \le 3.6),$$\end{document}(-0.9≤Q≤3.6), eckert number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(0 \le Ec \le 1.2),$$\end{document}(0≤Ec≤1.2), variable viscosity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${( - 100} \le \theta_{e} \le 100)$$\end{document}(-100≤θe≤100) and thermal conductivity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${(0} \le \varepsilon \le 4.1).$$\end{document}(0≤ε≤4.1). The repercussions of joule heating, wall suction, heat generation & absorption, magnetic field, viscous dissipation accompanying with variable characteristics of the fluid are also examined as well. Kinetics of viscous fluid with variable characteristics of fluid having solid body rotation over a permeable disk (having cylindrical geometry) are analyzed. We transformed the governing equations of heat transfer (accompanied by variable properties) and fluid motion in to self-similar non-dimensional differential equations by using the Von-Karman variables which are then further analyzed numerically by utilizing Adams Bashforth method. For a physical insight, results are manifested to scrutinize the behavior of velocity and temperature profiles for different emerging parameters graphically. Moreover, the values of nusselt number & skin friction co-efficient are also computed and physically explicated for the assorted parameters. Outcomes of current investigations are compared with prior work, to ensure the authenticity of the numerical method, and strong agreement is noted.

Three-dimensional flow problems produced by virtue of a rotating viscous fluid are experienced in a number of engineering-relevant, product-designing, industrial and scientific products and procedures for instance centrifuges, viscometers, jet engines, turbines, mixing in chemical chambers, vacuum cleaners, disk drives of computers and turbo-machinery etc. A classicistic problem in the field of fluid mechanics is the problem of fluid flow that is generated by a disk which possesses large diameter and is revolving about it's origin under the influences of constant angular velocity. Many years ago, Von-Karman 1 composed and resolved this problem. Anyone doesn't visualized an accurate circular flow close to the disk in this fluid problem since there is not applied pressure gradient that supplied the radial component of acceleration inwards and, as a consequence, fluid scrolled outside. Von-Karman 1 who was capable to transform the system of Navier-Stoke's equation into a set of ordinary differential equations with the assistance of appropriate similarity transformations. Cochran 2 noted some imprecisions in the solution of Von-Karman's problem and then developed an asymptotic solution of that problem. Ackroyd 3 elongated the problem of Von-Karman by accounting the impacts of injection/suction and illustrated series solution that comprises of the exponentially decomposing expressions. A similar problem is originated when the fluid revolves under the influences of constant angular velocity over a motionless plane is second fascinating problem. On the contrary, the fluid particles close to the disk moves in the direction of axis of rotation since it's circumferential constituent of velocity is decelerated under the similar pressure gradient. Tabassum and Mustafa 4 modeled numerically the heat transfer and partial slip flow of reiner-rivlin non-newtonian fluid generated by a revolving disk. A distinctive illustration made by researchers in this regard is given in references [5][6][7][8][9][10] . Contrary to Von-Karman 1 problems, Bodewadt 11 demonstrated that problem in 1940 by constructing boundary layer conditions. In Bodewadt's flow, fluid particles around the disk under the impact of radial pressure gradient directed to proceed inside towards the rotational axis where fluid is in the condition of stiff-body rotation. Bodewadt flow is monitored in a tea cup when rotation in the fluid is urged by continuously stirring and then left the flow of fluid for a little while. Bodewadt flow is additionally perceived further in many procedures such as hurricanes, foodprocessing industries, tornadoes, turbo-machinery and chemical mixing chambers etc. Rafiq et al. 12 examined numerical assessment of nanofluid for Bodewadt slip flow above a convectively heated penetrable disk. Sahoo et al. 13 obtained numerical solution of viscous fluid for steady Bodewadt flow using finite difference and keller-box method. Muhammad et al. 14 formulated entropy production in three-dimensional Bodewadt flow of nanofluid over an infinite stationary disk. Mukherjee and Sahoo 15 discussed the impact of radial stretch of the lower disk on the boundary-layer Bodewadt flow in the existence of the Coriolis force. Hina et al. 16 gives the numerical solution of rotating Bodewadt flow above a porous disk of micropolar fluid. Some recent attempts related to Bodewadt flow can be seen in Refs. [17][18][19][20][21][22] . Heat transfer is the field of thermal engineering which deals with the exchange, conversion, generation and use of thermal energy among physical systems (i.e.: heat energy). Transfer of heat is categorized within various mechanisms such as radiation, convection and conduction. Transfer of heat has broad range of applications in climate engineering, architecture, industries of chemical process, greenhouse effect and heat transfer in the body of human beings. There is an immense impact of temperature gradients on characteristics of fluid. Particularly, there is a direct relationship between temperature and viscosity in gases while an inverse relationship between temperature and viscosity is in liquids. Variable characteristics of fluid flows are often perceived in applications associated with aerodynamics & aero-acoustics. Takhar et al. 23 examined the impacts of fluid properties that depends upon temperature on boundary layer flow above a continuous moving surface. They considered two cases, namely, variable fluid viscosity and constant fluid properties. Salahuddin et al. 24 examined the rotating behavior of 2nd -grade fluid with mass & heat transfer effects between two parallel plates. Rafiq et al. 25 formulated the impacts of heat transport in fluid flow produced by a revolving permeable disk with variable characteristics of fluid. Ahmed et al. 26 numerically modeled the unsteady MHD flow and heat transfer with variable viscosity of nanofluid in carbon nanotubes over a shrinking porous surface. Adnan et al. 27 analyzed numerically the heat transfer in the nanofluid composed by silver nanomaterials and nanodimaond. Several attempts in this regime has been made by investigators includes [28][29][30] .
The heat which is produced due to the flow of an electric current through a conductor is called joule heating, also known as ohmic heating. Joule heating have a vast range of applications in our daily life like as glowing of filament of an incandescent light bulb, electrical fuses, hotplate of electrical tabletop etc. Viscous dissipation is an irreversible procedure in which the work is done by a fluid on adjoining layers owing to the action of shear forces is converted into heat. Repercussions of viscous dissipation plays a crucial role in natural convection in numerous devices that are subjected to big variations of gravitational force or which operate at high speed. Jawad et al. 31 discussed the analysis of transfer of heat and entropy generation in MHD flow of carbon nanotubes with the impacts of viscous dissipation and thermal radiation. Afridi et al. 32 35 gives the numerical solution of MHD flow of micro polar fluid under the impacts of joule heating and viscous dissipation over a shrinking sheet. Rasheed et al. 36 examined the MHD flow of chemically reactive casson liquid with effects of viscous dissipation and heat source over a penetrable stretching surface. The latest investigations in this regime include studies [37][38][39][40] . Fluids play dominant role in the application of flow rate and heat transfer in numerous engineering systems, biomedical procedures and technological developments. Gases and liquids are known as fluids because they can be made to move, or flow. The molecules of any fluid themselves are in constant, random movement, colliding with each other and with the walls of the container. A fluid is said to be viscous or real fluid if it doesn't flow easily because these fluids have finite values of viscosity that apply the shearing stress force on a surface with which is during a contact. These fluids have more resistance to flow. These fluids have very broad range of applications in our practical life. Equipment used for manufacturing requires suitable lubrication to run smoothly. The oil used as a lubricant for parts of heavy machinery should have a high viscous coefficient. Circulation of blood via veins and arteries rely on the viscosity of fluids. Few of the viscous fluids add texture to foods: for example, honey is quite viscous and can convert the "mouth feel" of a dish. The highly viscous liquid is utilized to damp the movement of some instruments and is used as brake oil in hydraulic brakes. Ayub et al. 41 analyzed numerically the free convection flow of MHD viscous fluid with heat generation and Newtonian heating over a rotating vertical plate. Yong-Min Li et al. 42 explained the mathematical modeling of entropy optimization due to a rotating cone in a convective viscous fluid flow.
Magnetohydrodynamics describe the study of movement of electrically conducting fluid (e.g. liquid metals and plasmas) in the presence of a magnetic field. The key hypothesis behind magnetohydrodynamics is that magnetic fields can generate current in a moving conductive fluid, which sequentially produce a force on the fluid and also alter the magnetic field itself. MHD have broad variety of applications in the field of Astrophysics, Engineering, Geophysics and Magnetic drug targeting. Shah et al. 43 analyzed the radiative effect radiation in magnetohydrodynamic casson nanofluid flow along with entropy generation and chemical reactions above a nonlinearily stretching sheet. Awais et al. 44 explained the mathematical modeling of prandtl MHD melted fluid flow towards an inclined cylindrical surface. Khan et al. 45 explained the entropy optimization rate in magnetohydrodynamic flow of viscous liquids with chemical reactions due to curved stretchable surface. Alreshidi et al. 46 discussed the thermophoresis effects and brownian movement in an incompressible and time independent flow of MHD nanofluid due to permeable revolving disk. Bibi et al. 47 presented theoretical analysis of magnetohydrodynamic carreau fluid above a revolving disk with the aid of von-karman transformations. Some remarkable applications regarding magnetohydrodynamics are cited here [48][49][50][51][52] .
By utilizing the knowledge of pre-mentioned literature, the objectives of present research is: • To find numerical solution of heat transfer effects in a Bodewadt flow in existence of joule heating, viscous dissipation, magnetic field and internal heat generation/ absorption. • Kinetics of viscous fluid with variable characteristics of fluid having solid body rotation above a permeable disk (having cylindrical geometry) are analyzed. • Utilizing the usual Von-Karman variables to convert the governing equations of heat transfer (with variable fluid properties) & fluid motion into self-similar differential equations which are then analyzed numerically by using Adams Bashforth method. • Results are manifested to scrutinize the behavior of velocity and temperature profiles for different emerging parameters with the assistance of graphs, tables & bar charts. • Outcomes of current investigations are compared with prior work, to ensure the authenticity of the numerical method, and strong agreement is noted.

Problem development and governing model
Let us assume the flow of viscous fluid with variable characteristics of fluid is in the condition of rigid body rotation over a penetrable disk through its origin as portrayed in the Fig. 1a, whereas the complete flow chart architecture is depicted via Fig. 1b. The disk rotates around vertical axis with uniform angular velocity ω . The components of velocity V r , V ϕ , V z are taken in the direction of cylindrical coordinates (r, φ, z) . The disk temperature is supposed constant on T w while the temperature of exterior flow is considered as T ∞ . A uniform magnetic field of strength B 0 is functional along the z-direction while the induced magnetic field is presumed to be negligible. The equation of motion, Navier-Stoke's equation and thermal energy equation with heat generation and absorption, ohmic heating & viscous dissipation takes the form as mentioned below: www.nature.com/scientificreports/ Here V r , V φ and V z in the above equations represents the constituents of fluid velocity in the radial, tangential and axial directions, (ρ) indicates the fluid density, C p indicates the specific heat capacity, (B 0 ) represents the magnetic field intensity, k(T) stands for thermal conductivity, (µ) indicates the dynamic viscosity, (T) denotes the temperature of fluid, (σ ) represents the electrical conductivity, (T ∞ ) represents the free stream temperature and (Q * ) stands for coefficient of heat generation/absorption . The stress tensor components are given as: We presume that viscosity (µ) manifests an inverse linear dependency of temperature.
Here (µ ∞ ) represents ambient viscosity of the fluid and γ denotes a constant. When we increase the temperature, the viscosity of the liquids decreases (γ < 0) whereas the viscosity of air and other gases increases (γ > 0) . Equation (7) may also be rewritten in the form as given below: We presumed thermal conductivity that is dependent on the temperature is of the shape (Tw−T∞) , in which K ∞ represents thermal conductivity at ambient. The suitable boundary conditions for the considered flow model are defined as: Note that the centrifugal force is well balanced by the radial pressure gradient in the exterior flow. Mathematically, we can write it as: Similarity transformation. Depending upon the scale of length v ∞ ω , we write the similarity variable

As the scale of velocity is
√ v ∞ ω and scale of time is ω −1 , therefore, we introduce the dimensionless constituents of velocity (F, G, H) & profile of temperature θ without dimensions as follows: By applying boundary layer approximations and then using formulas (10) in Eqs. (1-5); we acquire: www.nature.com/scientificreports/ The boundary conditions are converted into: In the above mentioned expressions (11)(12)(13)(14)(15), denotes the value of Prandtl number at the ambient, Ec = indicates the dimensional constant and Q = Q * ρc P ω represents the parameter of heat generation/absorption. Nusselt number and Skin friction component:. We define skin friction parameter as given below: where τ ϕ denotes the circumferential wall stress and τ r represents radial wall stress that can be evaluated by utilizing the Newtonian formulas as: By utilizing expressions (7 and 10) in expressions (16, 17 & 18), we obtain: Here Re r = rω 2 v ∞ indicates the rotational Reynolds number. Additionally, we define the local nusselt number as follows: represents the wall heat flux of the current research problem. By using formulas (10) Fig. 2c. However, the radial inward flow generated by the pressure gradient decelerates with an enhancement in the value of wall suction parameter A. The impact of suction leads to betterment in the frigidity rate of wall at a vast scale that is beneficial in numerous technical processes. The direction of the axial flow is directed upward in the absenteeism of wall suction, no matter how much greater the quantity of roughness parameter may be. The circumferential component of velocity G remains constant away from the disk while it is decreased significantly close to the disk due to the viscous drag as portrayed in Fig. 2b.
The axially independent radial pressure gradient compels the fluid particles close to the disk to move in the direction of rotational axis. For that reason, the direction of the axial flow in the Bodewadt flow is directed upwards as seen in Fig. 2c. The axially upward flow results into the non-physical solution of heat energy equation. Figure 2c illustrates that axial component of velocity becomes constant when we applied an appropriate value of suction velocity. All the three constituents of velocity manifest oscillations whose amplitude shrinks when axial distance increases. When A ≥ 1, the direction of the axial flow is inverted such that the fluid instead moves towards the solid surface. Depth of heat penetration also decreases when suction velocity increases. The results of Fig. 2c depicts that when larger A is accounted than greater quantity of (cold) fluid is transferred in the direction of the disk. As a result, layer of thermal boundary is suppressed with an increment in the velocity of wall suction and boosts transfer of heat from the disk as portrayed in Fig. 2d. Lastly, the axially upward flow is a consequence of an unphysical solution to the problem of advection-diffusion. To overcome this physical incompatibility, we can reverse the direction of the axial flow through the application of sufficiently high suction (15)  www.nature.com/scientificreports/ velocity. Therefore, it can be stated that the effect of damping arising as a result of wall suction can overcome the absolute or convective type instabilities in Bodewadt flow known from the published literature. Figure 3a, b includes the curves of temperature (θ) obtained by changing the value of prandtl number in two cases (A = 5 & A = 8). It is suggested from the comparison of Fig. 3a, b that when velocity of wall suction increases than the deepness of heat penetration is much suppressed. It can be noted instantly that an enhancement in the values of prandtl number rises the temperature of the fluid in the vicinity of the disk. The graphs demonstrate that the thickness of layer of thermal boundary and fluid temperature decreases with an increment in the prandtl number on a fixed value of ζ . The reason behind this is that prandtl number (Pr) is termed as fraction of momentum diffusion to thermal diffusion. When thermal diffusivity reduces, automatically the value of prandtl number (Pr) increases and fluid temperature decreases. In other words, temperature and prandtl number have an inverse  In heat transfer problems, the prandtl number controls the relative thickness of the thermal and momentum boundary layers. The prandtl numbers of gases are approximately equal to 1, which denotes that both heat and momentum dissipate through the fluid at about the same rate. Heat diffuses very slowly in oils (Pr > > 1) and very quickly in liquid metals (Pr < < 1) relative to momentum. As a result layer of thermal boundary is much thinner for oils and much thicker for liquid metals relative to velocity boundary layer. However, the quantity of θ ′ (0) increases considerably (in absolute sense) signaling a betterment in the rate of heat transition from the disk. The influences of magnetic interaction parameter (M) on the velocity constituents (F, G, H) in (r, φ, z) directions and temperature curves ( θ ) can be visualized in Fig. 4a-d. The radial, circumferential and axial components of velocity are decreasing gradually with intensification of strength of magnetic field (M) as evinced in Fig. 4a-c. The physics behind this is that an increase in the magnetic parameter (M) has the impact of damping the velocity profiles of fluid. This is because the direction of flow that is normal to the transverse magnetic field will result in a resistive force named as Lorentz force that is identical to the drag force which resist the fluid flow and as a consequence reduces the radial (F), axial (H) and tangential (G) velocities. Behavior of temperature curves (θ) for non-identical values of magnetic field parameter M is illustrated in Fig. 4d. Contrary to the profiles of velocity, temperature (θ) enhances when we increase the value of M. This is due to the fact that bigger drag co-efficient on the surface owing to sturdy magnetic field offers resistance to the particles of fluid. As a consequence produces heat owing to which temperature (θ) enhances. In modern technology, magnetic fields are especially used in electro mechanics & electrical engineering. Magnetic fields of rotating behavior are utilized in both electric generators and electric motors. Magnetic forces provide information about the charge carriers in a material via Hall effect. In electrical devices such as transformers, the magnetic field interaction is conceptualized and scrutinized as magnetic circuits. The Earth creates its own magnetic field, which protects the ozone layer of Earth from the solar wind and is significant in navigation using a compass. Figure 5a-d illustrates the change in velocity constituents (F, G, H) & temperature curves (θ) by altering the value of variable viscosity parameter ( θ e > 0). In these figures, the quantity of prandtl number Pr is kept persistent at Pr = 1 that is an adequate selection for air & numerous more gases. The variable viscosity parameter ( θ e ) elasticities the significance of the quantity 1 γ which is in direct proportion to the rate of dependence of viscosity on temperature of the fluid. Viscosity is a physical property that explores fluid's internal resistance towards its motion www.nature.com/scientificreports/ and oppose the shear deformation rate. A fluid having low viscosity easily flows because its molecular makeup results in very small amount of friction when it is in movement. A fluid with large viscosity resists motion because its molecular makeup gives it a lot of internal friction. We suppose the viscosity of fluid to vary as an exponential function of temperature in the dimensionless form, whereas the constant value of viscosity coefficient far away from the disk, is the variable viscosity parameter. When ( θ e → ∞ ), the fluid viscosity becomes equivalent to the ambient viscosity whereas the viscosity of fluid becomes a sturdy function of temperature when θ e → 0 . Figure 5c depicts that the axial velocity expedites in the downward direction as θ e increases. As a consequence, this effect causes a reduction in the severity of the radial flow as illustrated in Fig. 5a. The value of entrainment velocity H (∞) is also increased indicating a growth in the fluid volume that is sucked in the direction of the disk. The circumferential component of velocity G accelerated in the upward direction when the value of θ e increases as shown in Fig. 5b. As a result, the depth of heat penetration is (slightly) diminished and rate of heat transition is augmented when θ e is increased as portrayed in Fig. 5d. In Fig. 6a-d, we calculate velocity constituents and temperature curves for various values of variable viscosity parameter θ e when ( θ e < 0). In such situation, the viscosity of the fluid has an inverse relationship with temperature. Different to the outcomes of Fig. 5c, the axial component of velocity decreases when we increase the value of θ e that in results to the thickening of layer of thermal boundary as illustrated in Fig. 6d. Figure 6b reveals that the tangential velocity accelerates in the upward direction when the value of θ e increases.
The influences of heat generation/absorption (Q) and Eckert number (Ec) are demonstrated in Fig. 7a-b. In Fig. 7a, negative values of parameter (Q) are representing the heat absorption (Q < 0) and positive values are corresponding to heat generation (Q > 0). It is obvious from the Fig. 7a that an enhancement in the heat generation produces an increase in the thickness of thermal boundary layer whereas an increment in the heat absorption yields a decrease in the thickness of thermal boundary layer. An increment in the absorption of heat implies that the heat is withdrawn by some means from the system. That is the reason that why an enhancement in the parameter of heat absorption results a reduction in the thickness of layer of thermal boundary and an increment in the rate of heat transport. Contrary to this, an increment in the positive values manifest strong heat generation that enhances the fluid temperature in accordance with the definition of heat generation. Figure 7b portrays the repercussions of Eckert number on the temperature profile. An increase in the temperature is observed for greater values of Ec. An Eckert number (Ec) is utilize to calculate the loss of energy during configuration of flow. It tells us the relationship between enthalpy difference and kinetic energy of liquid particles, and is used to characterize heat transfer dissipation. For an augmentation in the values of Eckert number from 0 to 1.2, the Hence the temperature of the liquid increases. Figure 8a-d illustrates the change in velocity constituents and temperature curves for distinct values of parameter of variable thermal conductivity ( ε).Thermal conductivity is one of the thermo physical property. Thermal conductivity of a material is defined as the measure of its ability to conduct/transfer heat. Heat is transferred at a higher rate in materials of high thermal conductivity as compared to the materials which have low thermal conductivity. For example, metals are very effective at conducting heat because they have typically high thermal  www.nature.com/scientificreports/ conductivity, while the opposite is correct for insulating materials such as Styrofoam. That's why we use metal utensils in cooking to allow heat to move through it quickly and around what we are cooking. Figure 8c shows that volumetric flow rate and axial velocity are decreasing functions of thermal conductivity parameter ( ε ). Such decrease in the axial flow intensity is indemnified by an increment in the inward radial velocity near the disk as illustrated in Fig. 8a. A cross over in the curves of velocity constituent G shown in Fig. 8b revealing that azimuthal velocity (G) increases away from the disk and decreases near the disk when parameter ( ε ) increases. Figure 8d shows that temperature profiles (θ) becomes thicker when we change the values of parameter ( ε ). It is because of this reason that thermal conductivity is in direct proportion with the values of ε . Hence an increment in the values of parameter ( ε ) enhances the thermal condition of fluid that in result thickens the layer of thermal boundary. A comparative analysis has been carried out and results are displayed in Table. 1. A very good agreement is observed between existing results and those of Rafiq et al. 53 . Results of skin friction parameter and Nusselt number for distinct values of emerging parameters are presented in Tables 2 and 3. Further, bar charts displayed the tabular data of skin friction parameter and nusselt number for various parameters in Figs. 9, 10, 11, 12, 13, 14, 15, 16 and 17. Table 2 demonstrate the impact of A, ε, Q, Ec, M and θ e against skin friction parameter. In boundary layer flows, the coefficient of skin friction is a crucial non-dimensional parameter. It specifies the fraction of the  www.nature.com/scientificreports/ local dynamic pressure, 1 2 ρU 2 , that is felt as shear stress on the surface. Table 2 shows that value of skin friction parameter decreases when we increase the value of A and θ e whereas skin friction parameter increases when we increase the value of ε, Q, Ec and M . Table 3 highlighted the impact of A, ε, Q, Ec, M and θ e on (Nu). In fluid dynamics, the nusselt number is defined as the ratio of convective heat transfer to conductive heat transfer across a boundary in a fluid. Convection comprises of both diffusion (conduction) and advection (motion of fluid). When Nu = 1, it means that the heat transfer by pure conduction or we can say that fluid is stationary. The value of nusselt number between 1 & 10 is a feature of laminar flow or slug flow. The greater value of nusselt number corresponds to more active convection, with turbulent flow usually in the 100-1000 range. It is noticed that an increment in magnitude of A decays nusselt number while a direct trend is noticed for ε, Q, Ec, M and θ e .

Conclusions
In the current study a mathematical model has been presented for the influences of heat transport in a Bodewadt flow in the existence of viscous dissipation, wall suction, magnetic field, heat generation/absorption, joule heating and variable characteristics of fluid. The consequence of the current problem demonstrated that appropriate